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The perturbative QCD approach has achieved great success in the study of hadronic B decays. 
Utilizing the constrained parameters in these well measured decay channels, we study most 
of the possible charmless Bs — > PP, PV and VV decay channels in the perturbative QCD 
approach. In addition to the branching ratios and CP asymmetries, we also give predictions 
to the polarization fractions of the vector meson final states. The size of SU(3) breaking effect 
is also discussed. All of these predictions can be tested by the future LHCb experiment. 
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1 Introduction 

There is a. co ntinuous progress in the study of hadronic B decays since the so called naive factorization 
approachffl^ In recent years, the QCD factorization approach (QCDF)!^] and perturbative QCD factor- 
ization (pQCD) approach^l together with the soft-coUinear effective theoryl^l solved a lot of problems in 
the non-leptonic decays. Although most of the branching ratios measured by the B factory experiments 
can be explained by any of the theories, the direct CP asymmetries measured by the experiments are 
ever predicted with the right sign only by the pQCD approachP The LHCb experiment will soon run in 
the end of 2007. With a very large luminosity, it will accumulate a lot of Eg events. The progress in 
both theory and experiment encourages us to apply the pQCD approach to the charmless Eg decays in 
this wotW^ 

In the hadronic E{Es) decays, there are various energy scales involved. The factorization theorem 
allows us to calculate them separately. First, the physics from the electroweak scale down to b quark 
mass scale is described by the renormalization group running of the Wilson coefficients of effective four 
quark operators. Secondly, the hard scale from b quark mass scale to the factorization scale \/ Ams are 
calculated by the hard part calculation in the perturbative QCD approachP When doing the integration of 
the momentum fraction x of the light quark, end point singularity will appear in the coUinear factorization 
(QCDF and SCET) which breaks down the factorization theorem. In the pQCD approach, we do not 
neglect the transverse momentum kx of the light quarks in meson. Therefore the endpoint singularity 
disappears. The inclusion of transverse momentum will also give large double logarithms In^fc^ and In^a; 
in the hard part calculations. Using the renormalization group equation, we can resum them for all loops 
to the leading order resulting Sudakov factors. The Sudakov factors suppress the endpoint contributions 
to make the calculation consistent.'^ 

The physics below the factorization scale is non-perturbative in nature, which is described by the 
hadronic wave functions of mesons. They are not perturbatively calculable, but universal for all the decay 
processes. Since many of the hadronic and semi-leptonic E decays have been measured well in the two B 
factory experiments, the light wave functions are strictly constrained. Therefore, it is useful to use the 
same wave functions in our Eg decays determined from the hadronic E decays. The uncertainty of the 
hadronic wave functions will come mainly from the SU(3) breaking effect between the Eg wave function 



Table 1: The branching ratios and CP asymmetry calculated in PQCD approach, QCDF and SCET approaches 

together with Experimental Data. 
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-6.7 ± 16 


30 ±6 


39 ± 15 ± 8 
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and B wave function.' In practice, we use a little larger uji, parameter for the Bg meson than the B^ 
meson, which characterize the fact that the light s quark in Bs meson carries a littler larger momentum 
fraction that the d quark in the Bd meson. 

2 Results and Discussion 

For Bs meson decays with two light mesons in the final states, the light mesons obtain large momentum 
of 2.6GeV in the Bg meson rest frame. All the quarks inside the light mesons are therefore energetic and 
collinear like. Since the heavy b quark in Bg meson carry most of the energy of Bg meson, the light s 
quark in Bg meson is soft. In the usual emission diagram of Bg decays, this quark goes to the final state 
meson without electroweak interaction with other quarks, which is called a spectator quark. Therefore 
there must be a connecting hard gluon to make it from soft like to collinear like. The hard part of the 
interaction becomes six quark operator rather than four. The soft dynamics here is included in the meson 
wave functions. The decay amplitude is infrared safe and can be factorized as the following formalism: 



C{t) X H{t) X $(a;) x exp 



-s{P,b)-2 f ^7,(«.(Ai)) 



(1) 



where C(t) are the corresponding Wilson coefficients of four quark operators, $(a;) are the meson wave 
functions and the variable t denotes the largest energy scale of hard process H, which is the typical 
energy scale in PQCD approach and the Wilson coefficients are evolved to this scale. The exponential 
of S function is the so-called Sudakov form factor resulting from the resummation of double logarithms 
occurred in the QCD loop corrections, which can suppress the contribution from the non-perturbative 
region. Since logarithm corrections have been summed by renormalization group equations, the above 
factorization formula does not depend on the renormalization scale /i explicitly. 

The numerical results of the Bg decays branching ratios and CP asymmetry parameters are dis- 
played in reiP In all the decay channels for charmless Bg decays, only several are measured by the CDF 
collaborationP We show those channels together with results by QCDpl^ and SCET approaches!^ in 
table [TJ From those comparison, we notice that the measured branching ratios are still consistent with 
the theoretical calculations. Like the case in B decays, the calculated branching ratios from the three 
kinds of methods overlap with each other, considering the still large theoretical uncertainties. A global 
fit is useful when we have enough measured channels. 

In table [1] the only measured CP asymmetry in Bg K~Tr~^ decay prefer our pQCD approach 
rather than QCDF approach. This is similar with the situation in B decays. The direct CP asymmetry is 
proportional to the sine of the strong phase difference of two decay topologie The strong phase in our 
pQCD approach is mainly from the chirally enhanced space-like penguin diagram, while in the QCDF 
approach, the strong phase mainly comes from the virtual charm quark loop diagrams. The different 
origin of strong phases gives different sign to the direct CP asymmetry imply a fact that the dominant 
strong phase in the charmless decays should come from the annihilation diagrams. It should be noted 
that the SCET approach can not predict the direct CP asymmetry of B decays directly, since they need 
more experimental measurements as input. However, it also gives the right CP asymmetry for Bg decay 
if with the input of experimental CP asymmetries of B decays, which means good SU(3) symmetry here. 

For the Bg VV decays, we also give the polarization fractions in addition to the branching ratios 
and CP asymmetry parameters?! Similar to the B VV decay channels, we also have large transverse 
polarization fractions for the penguin dominant processes, such as Bg (pcf), Bg K*^K*^, K*^K*'^ 
decays, whose transverse polarization fraction can reach 40-50%. 



3 SU(3) breaking effect 



The SU(3) breaking effects comes mainiy from tire Bs{Bfi) meson decay constant and distribution am- 
plitude parameter, light meson decay constant and wave function difference, and various decay topology 
differences. As an example we mainly focus on the decays B — > tttt, B Ktt, Bg Ktt and Bg KK, 
as they can be related by SU(3)-symmetry. A question of considerable interest is the amount of SU(3)- 
breaking in various topologies (diagrams) contributing to these decays. For this purpose, we present 
in Table [2] the magnitude of the decay amplitudes (squared, in units of GeV^) involving the distinct 
topologies for the four decays modes. The first two decays in this table are related by U-spin symmetry 
{d — + s) (likewise the two decays in the lower half). We note that the assumption of U-spin symmetry 
for the (dominant) tree (T) and penguin (V) amplitudes in the emission diagrams is quite good, it is less 
so in the other topologies, including the contributions from the VF-exchange diagrams, denoted by £ for 
which there are non-zero contributions for the flavor-diagonal states tt'^'tt" and K^K^ only. The U-spin 
breakings are large in the electroweak penguin induced amplitudes Vsw, and in the penguin annihilation 
amplitudes Va relating the decays Bd — > K'^tt' and Bg K^K~ . In the SM, however, the amplitudes 
Vsw are negligibly small. 



Table 2: Contributions from the various topologies to the decay amplitudes (squared) for the four indicated 
decays. Here, T is the contribution from the color favored emission diagrams; V is the penguin contribution from 
the emission diagrams; £ is the contribution from the W-exchange diagrams; Va is the contribution from the 
penguin annihilation amplitudes; and Vew is the contribution from the electro-weak penguin induced amplitude. 
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In B^ ^ K 7r+ and if+7r , the branching ratios are very different from each other due to 

the differin g s trong and weak phases entering in the tree and penguin amplitudes. However, as shown 
by Gronau^i^, the two relevant products of the CKM matrix elements entering in the expressions for 
the direct CP asymmetries in these decays are equal, and, as stressed by Lipkinl^ subsequently, the 
final states in these decays are charge conjugates, and the strong interactions being charge-conjugation 
invariant, the direct CP asymmetry in K^n^ can be related to the well- measured CP asymmetry 

in the decay Bj] — > K~^t:~' using U-spin symmetry. 

Following the suggestions in the literature, we can define the following two parameters: 

„ ^ \A{Bg^n+K-)\'-\AiBs^7:-K+)\^ 
' - \A{Ba^7r-K+)\^-\AiBd-^7r+K-)\^' 

^ ^ Afi-p{Bd^n+K-) BRjBg^n+K-) rjBd) , 
A'^"p{Bg^TT-K+) BR{Bd^T:+K-) riBg)' 

The standard model predicts i?3 = — 1 and A = if we assume L/-spin symmetry. Since we have a 
detailed dynamical theory to study the SU(3) (and U-spin) symmetry violation, we can check in pQCD 
approach how good quantitatively this symmetry is in the ratios R3 and A. We get R3 = -OMtoH 
and A = —0.03 ±0.08. Thus, we find that these quantities are quite reliably calculable, as anticipated on 
theoretical grounds. SU(3) breaking and theoretical uncertainties are very small here, because most of the 
breaking effects and uncertainties are canceled due to the definition of R3 and A. On the experimental 
side, the results for R^ and A are:^ 

i?3 = -0.84 ±0.42 ±0.15, A ==0.04 ±0.11 ±0.08. (4) 

We conclude that SM is in good agreement with the data, as can also be seen in Fig. [T] where we plot 
theoretical predictions for i?3 vs. A and compare them with the current measurements of the same. The 
measurements of these quantities are rather imprecise at present, a situation which we hope will greatly 
improve at the LHCb experiment. 
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Figure 1: R3 vs A: The red (smaller) rectangle is the pQCD estimates worked out in this paper. The experimental 
results with their ±1(t errors are shown as the larger rectangle. 



4 Summary 

Based on the kx factorization, pQCD approach is infrared safe. Its predictions on the branching ratios 
and CP asymmetries of the B^{B^) decays are tested well by the B factory experiments. Using those 
tested parameters from these decays, we calculate a number of charmless decay channels —s- PP, 
PV and VV in the perturbative QCD approach. The experimental measurements of the three Bs decay 
channels are consistent with our numerical results. Especially the measured direct CP asymmetry of 
Bs — > 7r~K~^ agree with our calculations. We also discuss the SU(3) breaking effect in these decays, 
which is at least around 20-30%. We also show that the Gronau-Lipkin sum rule works quite well in the 
standard model, where the SU(3) breaking effects mainly cancel. 
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